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Abstract. The final stages of the evolution of electron-degenerate ONe cores, resulting from carbon burning in 
"heavy weight" intermediate-mass stars (8Mq ^ M ^ IIMq) and growing in mass, either from carbon burning 
in a shell or from accretion of matter in a close binary system, are examined in the light of their detailed chemical 
composition. In particular, we have modelled the evolution taking into account the abundances of the following 
minor nuclear species, which result from the previous evolutionary history: 12 C, 23 Na, 24 Mg, and 25 Mg. Both 
23 Na and 25 Mg give rise to Urea processes, which are found to be unimportant for the final outcome of the 
evolution. 24 Mg was formerly considered a major component of ONe cores (hence called ONeMg cores), but 
updated evolutionary calculations in this mass range have severely reduced its abundance. Nevertheless, we have 
parameterized it and we have found that the minimum amount of 24 Mg required to produce NeO burning at 
moderate densities is ~ 23%, a value exceedingly high in the light of recent evolutionary models. Finally, we have 
determined that models with relatively small abundances of unburnt carbon (X( 12 C)~ 0.015) could be a channel 
to explosion at low to moderate density (~ 1 x 10 9 g cm -3 ). This is clearly below the current estimate for the 
CZJ . explosion/collapse threshold and would have interesting consequences. 
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1. Introduction effects, such as rotation-induced mixing and overshoot- 
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their origin in the effects of rotation during the core hydro- 
gen and helium burning phases. Additionally, they explic- 
itly mention that the advanced stellar evolutionary stages 
appear too short in time to allow the rotational instabil- 
ities to have a significant impact during the late stages. 
Hence, given the lack of accurate multidimensional cal- 
culations for this phase of the evolution of heavy-weight 
intermediate mass stars, one should conclude, in princi- 
ple, that spherically-symmetric evolutionary calculations 
can still provide a reasonable approximation to the real 
behavior — see the recent review of Maeder & Meynet 
(2000) for a thorough discussion of this topic. 

These objects were long ago proposed as clear can- 
didates for the so-called accretion-induced collapse - 
a process first described by Canal & Schatzman (1976). 
Currently, it is widely accepted that ONe electron- 
degenerate objects do collapse to form a neutron star when 
the central density is driven beyond a certain threshold - 
see, for instance, Gutierrez et al. (1966, 1997), or Canal & 
Gutierrez (1997), and references therein. In particular, in 
the pioneering work of Miyaji et al. (1980) it was shown 
that all stars in the mass interval 8Mq ^ M <, 12Mq 
would develop electron-degenerate ONeMg cores during 
carbon shell burning. These cores would undergo a phase 
of electron capture on 24 Mg and 24 Na first, and later on 
20 Ne and 20 F, to finally ignite Ne and O explosively at 
central densities higher than ~2x 10 10 g cm -3 . At these 
very high central densities, fast electron captures occur- 
ing on the nuclear statistical equilibrium (NSE) material 
would rapidly drive the Chandrasekhar mass below the ac- 
tual mass of the degenerate core and, consequently, grav- 
itational collapse would ensue. Later, Woosley, Weaver 
& Taam (1980) and Nomoto (1984) reduced the above 
mass range to 8Mq M ^ 10 Mq. Stars more massive 
than ~ 10 Mq would ignite Ne in a series of non-explosive 
flashes and ultimately proceed through the O-burning and 
Si-burning stages in a way similar to more massive stars. 
Subsequently, Timmes & Woosley (1992, 1994) found that 
the density threshold for an ONe core to collapse to neu- 
tron star dimensions was around ~ 4 x 10 9 g cm -3 . On the 
contrary, if NSE was reached at densities below this crit- 
ical density, the degenerate object would be completely 
disrupted. Finally, Gutierrez et al. (1996) analyzed the 
role of Coulomb corrections, both in the equation of state 
and in the electron-capture threshold energies, and found 
that explosive NeO ignition takes place at densities high 
enough to ensure gravitational collapse to nuclear matter 
densities. The previously described scenario should apply 
not only to the cores of Super-Asymptotic Giant Branch 
(SAGB) stars (Garcfa-Berro & Iben 1994) but also to 
ONeMg white dwarfs accreting material from a compan- 
ion in a close binary system (Miyaji et al. 1980; Nomoto 
1984, 1987). 

Most of the studies done so far for this mass range 
had considered 24 Mg as a major constituent of the core. 
However, the detailed studies of Ritossa, Garcfa-Berro 
& Iben (1996), Garcfa-Berro, Ritossa & Iben (1997), 
Iben, Ritossa & Garcfa-Berro (1997) and Ritossa, Iben 



& Garcfa-Berro (1999) have clearly shown that the mass 
fraction of 24 Mg is much smaller than previously thought. 
Following the trends shown by these modern evolutionary 
calculations for the stars in this mass range, in this pa- 
per we investigate the effects of the chemical composition 
of the degenerate core on the possible outcome of these 
objects. In particular, we analyze the impact of the abun- 
dance of 24 Mg, the presence of 23 Na and 25 Mg, and the 
presence of unburnt 12 C, which appears to be a specific 
characteristic of current evolutionary models. The paper 
is organized as follows. In section 2 we describe the most 
relevant input physics and the initial model. Section 3 is 
devoted to the effects of an increasing 24 Mg abundance, 
while in §4 we analyze the role played by the previously 
unburnt 12 C. The role of the minor chemical species ( 23 Na 
and 25 Mg) is analyzed in §5. In §6 we discuss our results 
and we draw our conclusions. 



2. Input physics and initial model 

For the sake of conciseness, and given that the numeri- 
cal code and most of the physical inputs employed in the 
calculations discussed in this paper have already been de- 
scribed in previous papers — see Gutierrez et al. (1996) 
and references therein — we will not present them here in 
detail. Instead, in this section we only summarize the most 
relevant physical ingredients which are relevant for the 
discussion of the evolutionary sequences analyzed below. 
These are the electron capture rates, Coulomb corrections 
and the prescription adopted for convective transport. 

Electron capture rates have been taken from Oda et 
al. (1994), which provides the most up-to-date, and phys- 
ically consistent determinations. This point seems already 
fixed, as the rates for sd-shcll nuclei have remained al- 
most unchanged since the work of Takahara et al. (1989). 
However, the data of Oda et al. (1994) are much richer 
for 24 Mg, 20 Ne and other sd-shell nuclei relevant for this 
work, such as 23 Na and 25 Mg. 

Coulomb corrections to the electron capture threshold 
have been included as in Couch and Loumos (1974). The 
net effect on the Coulomb interaction energy is to lower it, 
hence reducing the Fermi energy and increasing the cap- 
ture density threshold. The exact expression to compute 
the energy shift is 



AE thl = ii{Z + 1) - n(Z) 



(1) 



where fi(Z) is the chemical potential of the nucleus of 
charge Z, given by 



^ Z) = - kB T (l)M°- 9 + Cl (|) 1/3 + C2 (l) 

+ o(i/r)J 



Z\V3 



do 
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where T = Z 5 / 3 e 2 /a c k B T = 2.275 x 10 5 Z 5 /3( p y e )V3/T is 
the Coulomb coupling parameter, o c is the inter-electronic 
distance, Z is the average charge, and the four constants 
adopt the following values: C\ = 0.2843, c 2 = —0.054, 
d = —9/16 and d\ = 0.460. The rest of the symbols 
have their usual meaning. All the electron captures, as 
well as the Urea processes, are considerably affected by 
the change in the threshold. In the case of the Urea pro- 
cesses, the circumstances are especially interesting since 
the electron capture threshold will be increased, while 
the threshold for the (3 decay will decrease. At the densi- 
ties and temperatures relevant for this study, the change 
in the threshold energy amounts to roughly ~ 300 keV, 
or to an electron-capture density threshold increase of 
~ 5 x 10 8 g cm~ 3 , thus favouring the collapse. 

Another very important issue is what convective cri- 
terion must be adopted. In the paper of Miyaji et al. 
(1980), convection was assumed to set in according to 
the Schwarzschild criterion: a convective core started to 
develop at the onset of the electron captures on 24 Mg. 
Convective heat transport contributed to keep the tem- 
perature below that of explosive Nc ignition (<~ 2 x 
10 9 K) along core contraction until a central density 
p iS 2 x 10 10 g cm~ 3 was reached. That was questioned by 
Mochkovitch (1984), who argued that the steep gradient 
of the electron mole number, Y e , produced by the same 
electron captures that create the superadiabatic tempera- 
ture gradient, should stabilize the fluid against convective 
motion and a semiconvective region would form. Miyaji 
& Nomoto (1987) tested this situation by assuming semi- 
convective mixing to be negligible. They found that ex- 
plosive NeO burning was ignited at a central density 
~ 9.5 x 10 9 g cm~ 3 : heating by the electron captures 
and cooling by thermal neutrinos was then purely local. 
Thus, convective criteria are essential in ascertaining the 
final outcome of the evolution of compact ONe objects. In 
this paper we adopt a rather pragmatic approach. That 
is, we first consider the classical Schwarzschild criterion 
V < V a d. However, since this stability criterion assumes 
complete chemical homogeneity and this is not the case 
when electron-captures start, we also use the Ledoux cri- 
terion: V < Vl, where: 



V L = V, 



ad 




(3) 



To account for the convective mixing, we have em- 
ployed a simple but robust numerical procedure. When 
convection sets in, the chemical composition is homo- 
geneized in the region given by 



min(A miXjk 

? ^conv,k X At) 



(4) 



where A m ; X: k and w CO nv,k are, respectively, the mixing 
length scale and the convective velocity of a given mass 
shell and At is the time-step. So, we are assuming a high 
convective efficiency in both chemical transport and en- 
tropy transport. 



Our initial model, used in all the calculations, is a 
1.375 Mq electron-degenerate object, whose main chemi- 
cal constituents are 16 and 20 Ne, with the distribution 
obtained by Ritossa et al. (1996). Its central density and 
temperature are, respectively, p c ~ 1.0 x 10 9 g cm~ 3 
and T c c± 2.0 x 10 8 K. The core is discretized in up to 
1000 mass shells by an adaptive algorithm, which en- 
sures good resolution both in chemical composition and 
in entropy throughout the structure. The ONe object is 
slowly contracting due to spherical accretion at a rate of 
<~ 10~ 6 Mq yr _1 . Its mass is consistent with the mass 
of ONe white dwarfs, and is close to the Chandrasckhar 
mass, making this object prone to experience sudden in- 
stability leading to either disruption or collapse. 

3. The effects of the abundance of 24 Mg 

Since the pioneering work of Miyaji et al. (1980), the astro- 
physical nuclear reaction rates have undergone significant 
changes. In particular, the amount of 24 Mg present in typi- 
cal ONe cores has been reduced by almost an order of mag- 
nitude (Ritossa et al. 1996). Currently, it seems that the 
most important rates are already reasonably well known 
but, nonetheless, it would be interesting to ascertain the 
minimum abundance of Mg required to cause an early 
explosion of the degenerate core. For this, we have param- 
eterized the abundance of 24 Mg while keeping constant the 
ratio O/Ne. While this is not truly consistent, it serves well 
the point relevant here: how much 24 Mg is the minimum to 
engage a explosion at about 4 x 10 9 g cm" 3 , which would 
possibly lead to the complete disruption of the degenerate 
object. With the currently accepted abundances of 24 Mg 
(~ 3% by mass), the maximum temperature attained dur- 
ing the electron captures is only 3.8 x 10 8 K (Gutierrez et 
al. 1996), very far below the minimum needed to initiate 
the explosive oxygen burning and to drive the matter to 
NSE. The goal of this section is to determine the min- 
imum abundance of 24 Mg required to produce a central 
temperature of about 1.5 x 10 9 K. 

In Fig. 1 several evolutionary sequences with increas- 
ing abundances (3%, 10%, 15% and 25%) of 24 Mg are 
shown. We have restricted our study to the Ledoux 
case because the preclusion of convection implies a 
stronger heating of the central region, due to the rel- 
ative inefficiency of the radiative-conductive transport 
of energy Evolutionary calculations performed with the 
Schwarzschild criterion (Fig. 2) confirm this point, even 
taking into account the mixing induced by the convec- 
tive transport. As it can be seen in Fig. 1, the higher the 
abundance of 24 Mg, the higher the central temperature 
achieved during the phase of electron capture on this iso- 
tope. 

Table 1 shows the net effect of increasing the abun- 
dance of 24 Mg, both for the maximum central temperature 
and for the central density at ignition. The first trend ob- 
viously appearing in Table 1 is that failed explosions lead 
to higher collapse densities — see Fig. 1. This is due to the 
fact that electron captures proceed in conditions of con- 
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9.6 9.8 10 10.2 
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Fig. 1. Evolutionary sequences of an ONe electron- 
degenerate object for increasing abundances of 24 Mg. The 
abundances of 24 Mg are, respectively, 3%, 10%, 15% and 
25% by mass. The model with X( 24 Mg)=0.15 undergoes 
an off-center explosion and the model with X( 24 Mg)=0.25 
produces a low density explosion, leading to the complete 
disruption of the degenerate core. 
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Fig. 2. Evolution of an ONe object with the Schwarzschild 
criterion (solid line). The dotted line corresponds to a case 
in which convective mixing is inhibited. In the case with 
mixing, the effects of a convective Urea shell are evident. 
The peak corresponds to the point where the nuclides of 
mass 24 able to capture electrons (namely 24 Mg and 24 Na) 
are exhausted throughout the convective core. 



(="Mg) 


Tmax,c 

(10 8 K) 


Pign.c 

(10 9 g cm" 3 ) 


0.00 




9.410 


0.01 


3.282 


9.412 


0.02 


3.561 


9.422 


0.03 


3.841 


9.431 


0.04 


3.957 


9.682 


0.06 


4.605 


9.847 


0.08 


5.261 


10.221 


0.10 


5.957 


10.409 


0.15 


7.484 


8.310* 


0.20 


9.118 


5.228* 


0.25 


20.000 


4.514 



Table 1. Ignition density and maximum central temper- 
ature reached during the captures on 24 Mg and 24 Na, for 
several initial abundances of 24 Mg. The cases whose cen- 
tral densities are marked with an asterisk correspond to 
off-centre explosions. In all the calculations, the Ledoux 
convective criterion has been adopted. 



stant pY c ~ n c , and hence higher amounts of 24 Mg imply 
larger reductions of Y e . Consequently, the collapse density 
increases from 9.41 x 10 9 g cm' 3 , when there is no 24 Mg 
in the degenerate core, to about 10.409 x 10 9 g cm" 3 , 
when the amount of 24 Mg is 10%. The minimum abun- 
dance of 24 Mg for which NSE develops during the electron 
capture phase is about 15%, and the explosion is clearly 
off-center. The reason is simple: during the phase of elec- 
tron capture, part of the energy generated is released on 



the mass shells above that in which the electron capture 
reactions are proceeding; thus, the temperature of these 
shells is constantly increasing. There is also a second ef- 
fect: a reduction in Y c increases the rate of contraction 
and, thus, the gravitational enery release, leading to a still 
higher temperature of the material. This can be readily 
seen in Fig. 1, in which the slope of the evolutionary track 
in the log(p c ) — log(T c ) plane becomes steeper with in- 
creasing abundances of 24 Mg. Eventually, the energy gen- 
eration in one of these preheated shells is able to drive 
the temperature to the point at which oxygen burning 
can be self-sustained, and, in fact, accelerated to produce 
the transition to NSE. For even larger amounts of 24 Mg 
the mass shell in which the off-center explosion ensues is 
located closer to the center. The explanation is, again, 
rather straightforward, as the preheating for increasing 
abundances of 24 Mg is more extreme and the change in 
temperature caused by the electron captures on 24 Mg is 
also correspondingly higher. 

A central explosion occurs only when the amount of 
24 Mg approaches ^25%, which is not credible given the 
current nuclear reaction rates relevant for the production 
of 24 Mg during the carbon burning phase. The same can 
be said about the minimum abundance needed to produce 
off-center explosions: it is extremely unlikely — and many 
things would need to change significantly, in particular 
the rates of nuclear reactions during carbon burning 
that 24 Mg could by itself induce the heating to NSE in a 
electron degenerate core at low density. 
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Fig. 3. Evolution on the log(p c )-log(T c ) plane for the 
models incorporating abundances of 12 C ranging from 
0.01 to 0.06. The only model that fails to produce an 
explosion at moderate densities is that corresponding to 
A( 12 C)=0.01. For all these calculations, the Ledoux cri- 
terion has been adopted. 
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Fig. 4. Same as figure 3, but using the Schwarzschild con- 
vective criterion. Only the model with X( 12 C)=0.01 fails 
to trigger an explosion at moderate densities. 

additional term to take it into account. The expression 
employed throughout this work is then 



4. Explosions induced by 12 C 

Recent models (Dominguez et al. 1993; Ritossa et al. 1996; 
Gil-Pons & Garcia-Berro 2001, 2002) have convincingly 
shown that carbon burning can be incomplete in degener- 
ate cores. The amount of remaining carbon depends some- 
what on the detail of the calculations but is typically of 
the order of ^ 1%. Nevertheless, Dominguez et al. (1993) 
found a central region of ~ 0.2 Mq with a high carbon 
content, of the order of 25%, that may eventually lead 
to an early explosion, and hence to a disruption of the 
electron-degenerate core. More specifically, Dominguez et 
al. (1993) followed the evolution of a 10 Mq star in a bi- 
nary system. The extreme dependence of neutrino cooling 
on density was able to cool down the center of the star 
to temperatures low enough to preclude the onset of cen- 
tral carbon burning. Although it has been claimed that 
this result could be to due the roughness of the numerical 
procedure adopted in Dominguez et al. (1993) — in par- 
ticular to the poor spatial resolution of their mesh — or to 
the very high mass-loss rate adopted, this is a point that 
undoubtely deserves more attention, since carbon burn- 
ing in degenerate conditions is prone to develop into the 
explosive regime, even when this element is found in only 
relatively small amounts. Woosley (1986) has estimated 
that, in the interval T 8 =2.5-7.5, for each 0.12% of carbon 
burned the temperature increases by 10 8 K. 

In this case, given that there is a difference on the 
atomic mass of 12 C and its daughter nuclei (namely, 16 O, 
20 Ne and 24 Mg), the Ledoux criterion must include an 



V L = V ad - 




'Ye 



(5) 



Given the strong dependence of the carbon burning 
reaction rates in density and temperature, it is expected 
that the ignition will ensue in a very limited region of the 
core. Hence, even when using the Ledoux criterion, convec- 
tion will develop at the very beginning of carbon burning. 
This will induce strong mixing on the timescale of the con- 
vective turnover. Hence, fresh material will replenish the 
carbon content of the burning region, effectively increasing 
the energy generation. Current evolutionary calculations 
predict a small amount of unburnt 12 C distributed more 
or less homogeneously within a substantial fraction of the 
core, with slightly higher abundances near the center and, 
consequently, the outcome of these calculations will de- 
pend sensitively on the adopted convective criterion. 

In order to determine the minimum amount of remain- 
ing carbon needed to completely disrupt the core at low 
densities, we have parameterized its abundance from 1% 
to 6%, a broad enough range. We have adopted both the 
Ledoux and the Schwarzschild criteria. Also, to investi- 
gate the worst possible case, we have artificially inhibited 
the convective motion and energy transport. In this case 
carbon burning proceeds in a strictly local mode, making 
the onset of NSE more difficult. This could be similar to 
what would happen if the thermonuclear timescale were 
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Convection 
criterion 


7m ax ( C) 

(10 s K) 


Pign(NcO) 

(10 9 g cm" 3 ) 


Outcome 


0.01 


L 


5.645 


9.671 


C 




S 


5.895 


18.643 


C 




_ 


4.242 


9.663 


C 


0.02 


L 


20.00 


0.537 


E 




S 


20.00 


0.537 


E 




_ 


6.901 


9.653 


C 


0.03 


L 


20.00 


0.450 


E 




S 


20.00 


0.448 


E 







9.084 


9.686 


C 


0.04 


L 


20.00 


0.365 


E 




S 


20.00 


0.363 


E 







11.22 


1.200* 


E 


0.05 


L 


20.00 


0.366 


E 




s 


20.00 


0.365 


E 






20.00 


1.165 


E 


0.06 


L 


20.00 


0.308 


E 




S 


20.00 


0.307 


E 






20.00 


1.108 


E 



Table 2. Models with variable 12 C abundance. In the table, "L" stands for the Ledoux criterion, while "S" stands 
for the Schwarzschild criterion; "C" denotes collapse and "E" disruption. For carbon abundances over 2%, all realistic 
models imply the complete disruption of the star. The case with a 4% abundance of carbon, and with the convection 
artificially inhibited (marked by an asterisk) corresponds to an off-center explosion. 



much shorter than the convective turnover and evolution- 
ary timescales. 

The first set of calculations is shown in Fig. 3, for which 
the Ledoux criterion was used. As it can be seen, the only 
model for which an early explosion is avoided is that with 
carbon abundance X( 12 C) ~ 0.01. In Fig. 4 the evolu- 
tion in the log(p c ) — log(T c ) plane is shown for the set 
of models in which the Swarzschild criterion was adopted. 
Again, only the models with X( 12 C) Ss 0.01 do not lead to 
an early explosion. Finally, in Fig. 5 the results obtained 
when convection was artificially inhibited are shown. In 
this case, those models in which the carbon abundance is 
<; 2% give an explosion at moderate densities. For com- 
pleteness, our results are summarized in Table 2, where 
the maximum temperature attained during carbon burn- 
ing (fourth column) , the ignition density (fith column) and 
the outcome of the evolution (sixth column) are shown as 
a function of the carbon abundance and the convective cri- 
terion. Table 2 shows that in the cases where convection is 
taken into account, the minimum abundance of 12 C should 
be as high as 1.5%, only a factor of 2 higher than the abun- 
dance presently found in most calculations. Even when the 
convective processes are artificially inhibited, the amount 
of 12 C required to achieve NSE at low densities increases 
to 4% — and the explosion occurs slightly off-center. It 
is clear that the presence of small amounts of unburncd 
12 C could be a possible channel for the complete disrup- 
tion of ONe degenerate objects at very low density unless 
non-physically consistent situations are artificially forced. 
It is also not surprising that a model adopting the 12 C dis- 
tribution obtained by Dommguez et al. (1993) results in 
an explosion at low density. Nevertheless, Rayleigh- Taylor 



instabilities would rapidly mix the 12 C-rich matter with 
the rest of the core. In this case, the abundance of 12 C 
is already within the range studied in the present work 
(about a 4.5%). 

Given the low density at which NSE proceeds, it is 
likely that some of the reactions will be rapidly quenched, 
thus producing an incomplete burning and a correspond- 
ing lower fraction of Fe-peak elements. While this is 
not compatible with classical type la SNe nucleosynthe- 
sis, it could be an alternative explanation for sublumi- 
nous type la SNe, like SN1991bg (Ruiz-Lapuente et al. 
1993). Clearly, this outcome deserves further considera- 
tion, which is beyond the scope of this paper. 

5. Effects of the Urea pairs 23 Na- 23 Ne and 

25 Mg .25 Na 

Recent evolutionary calculations including detailed nu- 
clear reaction networks have found that noticeable 
amounts of 23 Na and 25 Mg appear after the carbon burn- 
ing phase — see, for example, Gil-Pons & Garcia-Bcrro 
(2001). These nuclides have low electron-capture energy 
thresholds (Q =4.376 MeV and Q =3.833 MeV, re- 
spectively) and so they will start captures even before 
24 Mg, which has a higher energy threshold. The Urea pro- 
cesses for these nuclides are 23 Na(e,^) 23 Na(,e^) 23 Ne and 
25 Mg(e,^) 25 Na(,e^) 25 Mg, occurring at densities of 1.68 x 
10 9 g cm" 3 for 23 Na and 1.17 x 10 9 g cm" 3 for 25 Mg. The 
rates for these Urea pairs have also been taken from Oda 
et al. (1994). Note that throughout the calculations pre- 
sented in this section Coulomb corrections were applied to 
the rates, as was done previously. The chemical composi- 
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Fig. 5. Same as figures 3 and 4, but with convection ar- 
tificially inhibited. In this case, even a 2% abundance of 
12 C fails to trigger a moderate density explosion. 

tion of the ONe core was A( 16 O)=0.625, AT( 20 Ne)=0.262, 
A( 23 Na)=0.052, A( 24 Mg)=0.045, and X( 25 Mg)=0.015 as 
derived from Ritossa et al. (1996). Also note that for sim- 
plicity we have excluded the small abundance of unburnt 
12 C and renormalized the resulting abundances. 

As it can be seen in Fig. 6, the main effect of the Urea 
pairs consist of small decreases in T c and in Y c . The later 
effect is the only lasting one. These Urea pairs are never in 
equilibrium, as the core is slowly contracting and, hence, 
increasing the Fermi energy of the electrons. Some 10 9 
seconds after formation of the Urea pairs, the rates for 
electron-capture on the involved nuclides are much higher 
than the /?-decay rates. Eventually, all the 23 Na becomes 
23 Ne, and the 25 Mg becomes 25 Na, with a total reduction 
on Y c of 0.00286, which translates into a slightly higher 
ignition density for NeO burning. 

When using the Schwarzschild criterion, convection ap- 
pears during the electron captures on the nuclides of mass 
24 ( 24 Mg and 24 Na). Due to the limited extent of the con- 
vective region, there is no convective Urea for 23 Ne nor 
for 25 Na, since the upper layer of the convective region 
is well above the density for significant rates of [3 decay 
from these nuclides. Thus, the only durable effect of 23 Na 
and 25 Mg on the evolution of our objects is a small reduc- 
tion in Y c (from the original value of 0.4981 to 0.4952). 
Hence, due to this slight reduction in Y e , the ignition den- 
sity for NeO burning is 9.493 x 10 9 g cm -3 , clearly above 
the collapse threshold given by Timmes & Woosley (1992, 
1994). 

6. Conclusions 

The final fate of pure ONe electron degenerate objects 
seems at present well established: they will collapse to a 
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Fig. 6. Evolution in the log(p c ) — log(T c ) plane of an ONe 
electron-degenerate object with the abundances given by 
Ritossa et al. (1994). The effect of the Urea pairs involving 
23 Na and 25 Mg is apparent as a sudden decrease in T c near 
logp c =9.2. The convective criterion adopted was that of 
Ledoux. 

neutron star after electron captures on 20 Ne at a density of 
about 9 x 10 9 g cm -3 . The main goal of this work has been 
to elucidate whether the minor species present in ONe 
cores — namely, 12 C, 23 Na, 24 Mg and 25 Mg — can open 
a channel to produce an explosion at moderate densities 
(about 10 9 g cm~ 3 ) which would lead to the complete 
disruption of the star. 

We have found that the amount of 24 Mg required to 
produce a low density explosion is about one order of 
magnitude above the results of current, state-of-the-art 
evolutionary calculations. Hence, it is very unlikely that 
electron captures proceeding on 24 Mg (and its daughter 
nucleus 24 Na) can drive the temperature of the star to the 
point of igniting the NeO burning at p ~ 4.5 x 10 9 g cm~ 3 . 
Even in the (improbable) case that this actually happened, 
this density is near the current accepted threshold for 
gravitational collapse (Timmes & Woosley 1992, 1994). 
Consequently, if there is not a very significant change in 
the reaction rates leading to the production of 24 Mg, this 
nuclide can be discarded as a possible trigger of a moder- 
ate density explosion. 

On the contrary, even small amounts of unburned 12 C 
are enough to drive a moderate to low density explosion. 
The reason is twofold. First, the energy yield is much 
higher for carbon burning than for electron captures on 
the A = 24 nuclei. Secondly, there is much more 12 C avail- 
able, since the ignition is very local and can induce con- 
vective motion that leads to mixing that carries fresh 12 C 
into the burning region. In this way, even seemingly small 
abundances of 12 C (about a 1.5%) can produce an explo- 
sion at p c ~ 1 x 10 9 g cm~ 3 , which would lead to the com- 
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plete disruption of the star. This potential channel to com- 
plete disruption would give rise to a thermonuclear super- 
nova. The characteristics of such a kind of thermonuclear 
supernova are beyond the scope of this work, but could ex- 
plain in part why there seems to be a correlation between 
spiral arms and Type la supernovae (Bartsunov, Tsvetkov 
& Filimonova 1994; McMillan & Ciardullo 1996). Finally, 
the effects of the Urea pairs generated by 23 Na and 25 Mg 
are minor. Apart from a small reduction in the tempera- 
ture, the only lasting effect is a small decrease in the elec- 
tron mole number Y c . This will result in a slightly higher 
density for the electron captures on the nuclei of A = 24 
and A = 20. 

We note that multidimensional simulations of convec- 
tion indicate a tendency to increase the amount of convec- 
tive mixing and to induce rotational mixing to an extent 
that could slightly modify our conclusions. Many efforts 
are currently being invested in modeling carbon burning 
under strongly degenerate conditions using multidimen- 
sional codes — see, for instance, Reinecke, Niemeyer & 
Hillebrandt (2003) and Ropke, Hillebrandt & Niemeyer 
(2004a, b) for recent developments in this field. A time- 
dependent numerical procedure for mixing (necessary for 
coupling core contraction and convection) could slightly 
change the details of our conclusions, but it is unlikely 
to cause a qualitative change — that is, transforming a 
collapse into an explosion or vice versa. Until a complete 
treatment of multidimensional convection and burning is 
available, our results can be regarded as a fair approxima- 
tion to the real evolution of this type of star. 
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